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Abstract. Here we carry out computations that help clarify the La- 
grangian and Hamiltonian structure of compressible flow. The intent is 
to be pedagog ical and rigorous, providing concre te e xamples of the the- 
ory ou tlined in lHolm. Marsden. and Ratiul |l998(| and lMarsden. Ratiu. and WeinsteinI 
|1984|| . 
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1. Spaces and Actions 

We start with the group G = Diff(M) and the vector space V = C°°(M). 
MultipHcation in the group G is denoted by a single centered dot, i.e. given 
71,72 e G, 

71 ■ 72 = 71 072- 

For our purposes, we write the dual of V asV* = Den(M), the space of all 
C°° volume forms on M. If the manifold M has dimension n, then 

Den(M) = S1"(M), 

the space of smooth n-forms on M. We pair an element f €V with ji ^V* 
by 

(/Z,/)= / //i. 
JM 



1.1. G action on V . Let's define an action of G on 1/ by 

^^(/)=7*(/). 

What we have here is a map which associates a linear isomorphism <^Y^ to 
each element 7 G G. We will sometimes write this action more compactly 
using concatenation: 

7/ = $^(/) = 7*(/). 
The action is clearly a right action: 

(71 -72)/ = (71 o72)*(/) 
= 72*(7i*(/)) 
= 72(71/) 



1.2. Induced G Action on V* . The G action on V induces a right ac- 
tion on V* via the inverse of the dual of ^Y, . We now calculate the dual 
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isomorphism (^>^)*. Given f €z V and cu gV*, 

(($:^)*H,/) = (u;, 

= (^,7*(/)) 

JM 

= [ 7*(/7*H) 

JM 



= / /7*M 

JM 

= (7*(^^),/)- 



This shows that 



($^)*H=7*(a;) 
Now we can define the action of G on V*: 

= (7*)-'H 

= 7*H 

The proof that this is a right action of G on V* is identical to the case of G 
acting on V. 

1.3. Semidirect Product Group and Algebra. We will be concerned 
with the semidirect product space G®V* , which is a group consisting of the 
set G xV* equipped with the product 

■ (72, '^2) = (71 •72,^^2 + 72^1)- 

The Lie algebra g of G consists of the space TeG equipped with the ap- 
propriate bracket. As can be shown without much difficulty, for a given 
diffeomorphism 7 G G, 

T-yG = { all maps X : M ^ TM such that ttm o A = 7}. 

Therefore, 

TeG = X{M) = { all vector fields A : M ^ TM}, 

since ttm o A = idjv^ is precisely the condition that makes A a vector field. 
(Clearly e G G is the identity map idM on M.) 

Proposition 1. The bracket [^,?7]g on the Lie algebra Q is minus the Jacobi- 
Lie bracket of vector fields on M. In symbols, 
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Proof: Let denote the right-iiwaii&nt vector field associated with C,, 
i.e. X^i^j) = C ° 7- Keep in mind that X^ is a vector field on Diff(M). 

Let's verify that X^ is right-invariant. Set R-y{P) = (3°^, right-translation 
on Diff(M). Then, for all maps Ca ■ M ^ TM such that vrjv/ o = a, 

TaR-y • Ca = C o 7, 

It is now clear that 

T,R^ ■ X^irj) = (C o r/) o 7 

= (o7?o7 = Co R^{r]) 
= Xf oi?^(,?), 

as required. With this in mind, we use the local formula for the bracket 
(again, the Jacobi-Lie bracket of vector fields on Diff(M)): 

[X^,X^]{e) = DX^ie) ■ X^{e) - DX^{e) ■ (e) 

= DiO ■ 7] - D{r]) ■ C 

and this is the local formula for the Jacobi-Lie bracket of vector fields on 
M. We have therefore shown that 

(1) [X^,X^]ie) = [C,v]M. 

Now we are almost done. Typically, we define a bracket on q via 

(2) [C,vU:=[X^,Xf;]{e) 

where X^ is the Ze/t-invariant vector field associated with ^ G g. 
Lemma 1. Defining a bracket on the Lie algebra g as in ^ implies that 

(3) [Xf,X«](e) = -[e,r?],. 
Proof: Introduce two spaces: 

= Lie algebra of all right-invariant vector fields on Diff(M) 

= Lie algebra of all left-invariant vector fields on Diff(M) 

Start with the diffeomorphism on G given by (p^g) = g^^. One can show 
that for X a left-invariant vector field, f{X) = is a right-invariant 

vector field. Indeed the map / defined in this way is a Lie algebra isomor- 
phism from to X^. Furthermore one can compute 

for all E g. Now we have, for G g, 

f{X};){e) = {T^oXf;o^-^){e) 
= Tecl)oX^{e) 
= Te4> o (T(idG) o rf) 
= Te4> ■ r] 
= -V- 
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Then 

[X^,X^]ie) = [^,rj]M 

= /(^hC,-.l«)(-) 
= -K,-??]g 

which is what was required. Putting together (jSJ with finishes the proof 
of the proposition. 

1.4. The dual algebra g*. We wih take q* to be the space of one- form 
densities: 

Q* = {uj (g) fi : to e A^(M),/i e Den(M)} 
Then for,5G£|,^(g)/zG0* the pairing between the two elements is given by 



Jm 

1.5. Induced g action on V. The action of G on 1/ induces an action of 
g on V. We use <I> : G — > GL{V,V) to denote the original action: <I>(7) = 
7*( • ). Then we can compute Tg^ : Q — > L{V, V) as follows. Fix ^ € g and 
write ^ as a tangent vector to a curve 7*^ on G = Diff(M): 



^ de 



Y where 7*^ = id^. 



e=0 



Then for f eV, 

(T.*. ?)(/)= I 



(7-)"(/) = £(/. 

e=0 



by the dynamic definition of the Lie derivative. From now on we wiU use 
concatenation to denote the action of g on V, so that for ^ € g and f (zV, 

U = £d- 



1.6. Induced g action on V*. The g-action on V* is defined as minus the 
dual map (i.e. linear algebraic adjoint) of the g-action on V. We can also 
start with the action of G onV*; this action induces a g-action on V*. Both 
g-actions on V* are in fact the same, as we will now show. Take fj, G V*, 
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f eV, and 4 e 0. Then 

{fi,£^f)= f 
Jm 

= [ {£^{ffi) - f£^ix) 

JM 

= [ idi^if^i) + kdif,,))- [ f£^f, 

Jm Jm 

Since //i is an n-form, cZ(//x) = 0. Also, since dM = 0, Stokes' theorem 
imphes that 

/ die(//x) = 0. 
Jm 

Therefore, we have 

{fi,£^f) = {-£^n,f), 

so the 0-action on V* is 

^/x = £^ii. 

Verifying the calculation using the alternate approach, we start with the 
G-action on V* denoted by ^ : G GL(y*,V*), where for each j E G, 
= 7*(-). Then we compute Te^' : g L{V*,V*). As before, we fix 
^ G g and write 



de 



7^ where 7° = idc. 



e=0 



Then 



(7')'(».) = £(H 

e=0 



by the dynamic definition of the Lie derivative. 

1.7. Adjoint action (group). We start with the conjugation map 

I-yiv) =7-?7-7"^ 
We can then write the adjoint action as 

Ad^C = Tel^ ■ C 

= r7 o e o 7-1 

= 7*C 
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1.8. Dual of adjoint action (group). Next on our list is the Ad* action. 
We take a; (g) ;u G g*, ^ G g, and 7 G G. Then the calculation runs as follows: 

(Ad*a; 0n,O = A*, Ad^O 
= (w (2> /X, 7*0 

i^^jw IJ, 

M 

Jm 
Jm 

= / in ^7 M 

= (7*a;0 7>,O 
= (7*(^®m),0, 

Ad*(a; (8)/Lt) = 7*(a; (2>/x). 



so 



1.9. Adjoint action (algebra). Moving to the Lie algebra, we have an- 
other pair of adjoint actions. The adjoint action on the Lie algebra is simply 
the bracket 

ad^ : ^ 0,ad^r? = [^,r]]g = -[C,ri]M- 

1.10. Dual of adjoint action (algebra). Fix oj^fi & g* and ?? G fl. Then 

the calculation runs as follows: 

(ad|(u; ® ii),ri) = {lo ® fi, ad^r]) 



/ 



IM 

Now we apply the following identity: 

We apply Cartan's magic formula to the first Lie derivative: 

£^\rfOJ = di^i^o; + i^di^w 
Note that i^w is a function, so i^ijjo; = 0. Then the divergence theorem yields 

/ i^di^w = - / i,,^; div^^ ji. 
JM JM 

-/ ■[^,v]m IJ'= i^(£^a; + c<;diV/,0/x, 

JM JM 

ad|(a; (8> yu) = {£^i^ + i^drv/^^,) (8> /j. 



Hence 



IM Jm 

which means that 
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1.11. Diamond map. We start with the map pj. Given ^ G 5 and / G V^, 
we define 



1.12. Summary of Results. Starting with 

• 7 G G = DifF(M), 

• / G y = e°°(M), 

• peV* = Dcn(M), 

• ^) £ = ^{M), and 

• e®v eg* = Ai(M) ® Den(M), 

we summarize our work so far in the following table: 

G action on F : 7/ = 7*(/) 
G action on F* : 7// = 7* (//) 
Lie algebra bracket on g : r/Jg = —[^,rf\M 
Q action onV : = £^f 
Q action on V* : ^fj, = £^fj, 
Ad action of G on g : Ad-y^ = 7*^ 
Ad* action of G on g* : Ad* (6' (?) z/) = 7* (6* (?) v) 
ad action of g on g : ad^?7 = r]]g = — t^Jm 
ad* action of g on g* : ad|(^ ® u) = {£^9 + ^div,^^) (2> u 
xV* ^ Q* : f0H = df ^ H 

2. Unreduced Lagrangian Dynamics 




{p}{a),0 = {a,Pf{0) 
= {a, £^f) 




= {df0a,O 



Hence p*f{a) = d/ a. Now we can define the diamond map, 

o:V xV* ^g* 
/ o a = d/ (g) a 



Our goal now is to present a derivation of the Euler equations for com- 
pressible flow in material coordinates. First we write down an expression 
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for the unreduced Lagrangian. Here rj G Diff(M), fj € Tj^T)iS(M), and 
H^eV* = Den(M). Then the Lagrangian L : rDifF(M) x ^ M is 



M 



\Mm)\\'' -w[[r,,l,oMm))) 



Mo- 



The first term in the integrand is the kinetic energy, while the second term 
is the potential energy. We assume the fluid is barotropic, meaning that the 
pressure is a function only of the fluid's density. This forces the potential 
energy to also be a function only of the fluid's density. The density fiQ here 
is the material fluid density. By this we mean that //q is the density as a 
function of material points in the fluid, at time t = 0. 

2.1. Equations of motion. We obtain equations of motion by applying 
Hamilton's principle directly to L. The calculations that follow are almost 
identical to the calculations that would result from taking L and plugging 
it into the Euler-Lagrange equations. 

First we write mo = Po d"a; for some po G 6°^ (M) . Then we have 

dL= [ [ffSrt 
Jm 

- W ((det Vr/)-VoM) (-(det Vr?)-^(V77)-^ • VSrjpo)] pod^x 

[ffSrf 



M 



+ W ((detVr7)-Vo(m)) ( (det Vr/)-i^i£j^^po 



Po d"a;, 



and after integrating by parts, 

SL= [ {-Wpo 
Jm 

d 



dXj 



W' ((det V?7)-VoM) (detV?7) 



pI 



5r,' \ d"x. 



Along solution curves, 5L must be zero for all 5rf^ implying that 77 must 
satisfy the PDE 



(4) Poff 



d 



dX3 



W^'((detVr?)-VoM) ("(detVry)-^^ 



7-1 V 



■pI 



2.2. Substitution of variables. We will now try to express this PDE in 
the spatial variables u and ps- This will enable us to compare it to standard 

compressible fluid PDE. The spatial density ps (a function of spatial points 
X and time t) is related to the material density po in the following way: 

1 



177 



Differentiating the relation 



det V?7 
uo rj = rj 



po- 
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gives us 

7] = u o T] + Vu ■ uorj. 
Substituting into Q), we have 

[(det S/r])ps o r]]{u o 7] + Vu ■ uo7]Y 
We can now get rid of or] everywhere to obtain 



W' {ps o r/(m)) ipsOT] —. — {detVr])ps o t] 



(5) {detVr])ps{u + Vu- u) 



We have used the following fact: because = r]^{X), we have, locally, 

d dx^ d drj^ d 

We expand the right-hand side of © into three terms: 



The last term is actually zero: 
dr]'' d fdi-q-^y , , 1 



dXi dx^ \ dx^ ' det(V77-i) 



\dX^ dx^dx^ det(Vr/-i)y (det(Vr/-i))2 

~ \dXi dx^dx^ det(Vr?-i)y ~ (det(V??-i))2 9X™ ^ 
= 0. 

This means we can simplify ((SJ all the way to the following: 

+ Vn • uY = -^[W'{ps)ps] - W'iPs)^. 

2.3. Induced action of G on TG. Before we show that L is invariant 
under the action of the group G = Diff(M), we revisit the action of G 
on TG. For each 7 € G, let $g : G ^ G denote the action of 7 on G. 
Specifically, 

$^(ry) =7707. 
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This action induces an action on each tangent space T^G. Specificahy, we 
have the map r$-y : TG — > TG which acts as follows. Fix -q G T^G. We can 
write this element of T^G as = ^ I^^q where rf = rj. Then 



= rj oj. 

Hence we obtain the action of 7 on TG: 

7 • (77, rj) := • (77, 7)) = (ry o 7, 77 o 7). 

Or since 77 o 7 G r^Diff(M), we can ignore the effect of the action on the 
base point of the tangent vector and write 



The group G acts on Den(M) via pullback, i.e. 7/i = 7*/x. (Though some 
actions considered are actually right actions, we will always write the group 
variable on the left.) 

2.4. Proof of G-invariance of L. We now verify that our Lagrangian L 
is in fact G-invariant. 



r^$^ • 77 = - $^(77^) 



e=0 




7 • 77 = 77 O 7. 



L(77 7,7» 




Recall that 



(77 o 7)*(7» = ri^-f*j*n = ri^n, 



so that 




Noting that 7*(/) = / o 7, we see that 




Now applying the change of variables theorem, we have 




which is precisely L{fj,fj,), proving G-invariance. 
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3. Reduced Lagrangian Dynamics 



Now we can write down the reduced Lagrangian. Here v G X(M), and 
fis e Den(M). 

l{v,Hs) = L{vor],rfiis)- 
Note that the Eulerian velocity v = fjo r]~^, and w is a tangent vector based 
at the identity element of Diff(M). Also note that jig = ??*Mo or fj,Q = rf Hs- 
We write the right-hand side as follows: 



L{v or],r]*fj,s 



I f 

Jm . 



1,1 



Using the properties of pullback, we write 



V 



M 



Now by the change of variables theorem, 



l{v,iis) 



M 



1,1 , 

-\\v\ 
2" ' 



Wills) 



which is our final expression for the reduced Lagrangian i : g x y* ^ M 
where, of course, q x V* = X{M) x Den(M). In what follows, we will use 
the following notation: fig = Ps d"aj, where d^x is the canonical n-form on 

M and p., G e°°(M). 

3.1. Equations of Motion I: Variational Principle. Our first path from 
the reduced Lagrangian I to the equations of motion will involve the reduced 
variational principle. That is, critical points [v^Ps) of Z, subject to the 
constraints 

. dw . , , 
ov = -TT- + \v,w\, and 



Sps 



dt 



are actual trajectories of the system. Or, in other words, if one solves the 
variational equation 

Sl{v,ps) ■ {5v,dps) = 

subject to the above constraints on the variations Sv and Sps, one finds so- 
lutions to the equations of motion that describe the evolution of the system. 

3.2. Derivation of Constraints. Before proceeding, we confirm that the 
above constraints are natural. They arise precisely from the relationship 
between the unreduced (rj, po) and reduced (v, ps) variables. Free variations 
of the former correspond to constrained variations of the latter. 

To sec this, we start with the relationship between the reduced and unre- 
duced velocity variables, u = rj o r]~^, and calculate variations. First we 
write 



u 



rf o 
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5u := -f- 

de 



e=0 



f c- -1 T-, • d 

u = 01) orj + vrj ■ — 

de 



e\-l 



e=0 



Now because (77"^) ^ o (ry*^) = id, we can differentiate both sides with respect 
to e at e = 0. This yields the following formula: 



Te 



e=0 



Applying this to our earlier expression for 5u, we have 

Su = dr] o r]~^ — Vr) • Vr/~^ • Sr) o r]~^. 
Defining w = 6rj o r)~^ , we find that 

w = 61) o Tj"^ + V5r] ■ Triv^^)- 
ot 

We can calculate the very last term in this formula by differentiating rj^^or] = 
id with respect to t, resulting in 



dt 



rj ^ = —Vrj ^ ■ rj orj ^. 



Substituting into our earlier expression for w, wc have 

w = Srj o r] '^ — V5r] ■ Vr]~^ ■ f] o ri~^ 
= Srj o r]~^ — Vw ■ u. 

Using what we have so far, along with the fact that 
we have 

Su = w + Vw ■ u — Vn ■ w = w + [ujw]. 
Next, using basic facts about the pullback of volume forms, we write 

1 



Taking variations, 

dpsOri + Vps ■ Sri = 



det Vr? 



Po- 



Po 
det Vrj 



(Vr/)-'^ • VSr], 



where the final • denotes pairing of linear transformations. That is, A - B = 
tr(A-^B). Now we solve for Sps, obtaining 



Sps 



-Vps ■ w — /9str(Vr7 ^ ■ VSrj). 
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In coordinates, 



dxo dX 
d{Sr] o r)~^y 



Therefore, 



Sps 



-V ps ■ w — /9sdiv(ti;) 



= div(u)). 



-d\v{psw). 



3.3. Critical Points of I. We now vary {v,ps) subject to the constraints 
which have just been verified. 

^1, 



51 ■ {Sv,Sps) 



M 



[v'5v' - W'{ps)?>Ps\ Ps + 



\vV -w 



5ps \ dx 



Using the constraints and switching to coordinates. 



51 ■ {5v, 6ps 



I 

Jm 



+ 



dt dx^ 



d 



dxi 



d 



+ W'{Ps)q-{Psw')\Ps 



7,\A-W\—{psW^)\dx 



Integrating by parts, 
51 



{5v,5ps) = ^-^^{v'ps)w' - 



d ■ ■ ■ dv^ 
dxJ^ dxr 



1, , 



W{ps 



PsW^ > dx. 



The above must hold for all w, so we have the following equation (in what 
follows, we write p for ps and use subscripts to denote spatial derivatives): 



dv 



idp 



dp 



--— /9 — V^— V^jV' p — V^V-'-p — V^V^ — — r 



dt 



dt 



,P - [W'{p)pU + [v^vi, - W\p)p^,]p = 

The third term cancels with the sixth term. Also, the second, fourth, and 
fifth terms summed together equal the left-hand side of the continuity equa- 
tion, i.e. 



dp 



dp 



-7^ + iTaP + v'-^T—- 



0. 



dt dx^ 
After eliminating these terms, we cancel a factor of p and obtain 

dv^ 



dt 



[W'{p)p,^ + iW'{p)p), 



Not surprisingly, this is precisely the same equation we got by taking free 
variations of the unreduced Lagrangian and then changing from material to 
spatial coordinates. 
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3.4. Equations of Motion II: Euler-Poincare Equation. We can also 
derive the equations of motion from the reduced Lagrangian by making use 
of the Euler-Poincare equation from semidirect product reduction theory: 

— —--ad*— — o 
dt Sv ^ Sv Sfig 



3.5. Derivatives of I. In order to apply the Euler-Poincare equations of 
motion, we must first calculate the variational derivatives of 1{v,ijls)- First 



d 
de 
d_ 
de 



e=0 



l{ve,lJ-s) where ^ 



e=0 



JM ^ 

_ /• 1 _d 
Jm 2 de 

= [ {v,0\ 
Jm 



e=0 
■ WtJ-s 



where the b operation is pointwise with respect to (•, ■ 

v\0 = {v,Om- 



I.e. 



Hence 



SI 
Sv 



Note that W only depends on the ps "part" of /x^, not on the canonical 
volume form d"x part. Therefore, 

W{fis)=W{ps<Px) = W{ps). 

Here W : Den(M) Q°°(M). Using W, we can write an alternate expres- 
sion for BW : Den(M) ^ L(Den(M), e°°(M)) as follows. Let p = pdJ'x 
and p' = p' (F'x. Then 



T>W{p)-p' = 



d_ 
lie 

d_ 
de 

dW 



t=0 



W{p + ep') 

I 

Wip + ep') 



dp 



€=0 

ip)-p' = W'{p)p'. 
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Furthermore, for densities fi, ni, and fi2, we have 

{■DWifi)-fii,fi2) = (w'{p)pi,fi2 



M 



I 

Jm 



W'{p)pil22 

w'{p)pip2(rx 

W'{p)p2fll 



W'{p)p2,ni 
{-DWip)-p2,Pi) 



Armed with these facts, we proceed: 



/ SI \ _ ± 
\6ps' I de 



1{v,pI) where 



d 



e=0 



1,1 , 



de 

2 - W{p[ 



11% = a 



e=0 



\v\\^ -W{p) 



a + ps [-TyW{ps) ■ a] 



e=0 
d 

M de 
"1 
M [2 

1,, ,, 



-\\v\\''-W{ps)-^W{ps)-Ps,a 



W{ps),a)-{'DW{ps)-a,ps) 



so we must have 



_SJ_ 

5ps 



2IM 



W{ps) - W'{ps) ■ Ps. 



3.6. Plugging into Euler-Poincare. In what follows, we write p to mean 
Ps. In an earlier calculation, we established that for ^ G g and 9 E Q*, 

ad|(0 z/) = {£^9 + ^div^O ® ^■ 
Applying this in our case yields 

d SI 
— (t!^ <Si p) = -{£vv' + v'diYnv) ® p + — o p. 
at op 

Now because div^t; is a scalar function, the two objects —v^dw^v p and 
—v' ® divuvp are equivalent. Thus we can pick off the advection equation 



which simplifies to 



dp 
'dt 



-{divnv)p = -£^p. 
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Then we can simplify our earlier equation: 

^ fJ. = -IvV (81 // + d ( — I (g) /i. 

We have made use of our previous computation for the diamond map. Now 
we apply Cartan's magic formula and substitute for the remaining derivative: 

dv" 1 

— + d\J + \M = -d\J - dW{^i) - d{TiW{^i) ■ 11). 
A coordinate computation verifies that, in fact, 

]^d\,v^ + \M = [{v-V)vt 
Taking sharps of both sides, we end up with 

_^ + . V)t; = - [dWiii) + d{T>W{^x) ■ ^)]« , 
which agrees with our previous results perfectly. 

4. Hamiltonian Dynamics 
We will attempt a Hamiltonian derivation of the following system: 



dtp + (pu)^ = 

dt {pu) + {pu^ +p)^ = 0. 

These are the Euler equations for a compressible fluid on the line. The 
equations are also known as the continuity and momentum equations, re- 
spectively, of one-dimensional gas dynamics. 

4.1. Setup. We know we can write the one-dimensional Lagrangian as a 
function of a velocity u and a density p = pdx: 

(7) l{u,p) = ^ j p\\ufdx + j pW{p)dx 

A simple computation reveals that 51 /5u = pu, which is therefore the correct 
momentum variable to use in the Hamiltonian formulation. Let us set m = 
pu and write H as a function of m and p: 

(8) H{m,p) = ^ j ^\\mf dx + j pW{p)dx. 

Here H is a function on the semidirect product Lie algebra QxV. Typically 
we view as a function on the dual of a Lie algebra, but in this case it does 
not matter. The reason is that we can identify m with the one-form density 

in = pdx 

and we can also identify p with pdx. With these identifications, we can 
write H (in, pdx) and its expression will be identical to Now we can 
apply the semidirect product reduction theorem, obtaining the Lie-Poisson 
equations on the dual algebra (g x V)*. (Note to reader: I will make this 
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more concrete and explicit in the next draft. For now I will focus on the 
computations.) 

4.2. Brackets. The semidirect product Poisson bracket is given in our case 
by 
(9) 

dx ' ' " 



{F,G}(m,p) 



m 



5F_ 6G_ 



5G 



£5G/Sm-^ I dx. 



The bracket under the first integral sign is the Lie algebra bracket on q. It is 
a standard result that this bracket is minus the ordinary Jacobi-Lie bracket 
of vector fields on M. Therefore in one dimension we have 

(10) {F,G}(m,p) = / m — — -— — 



(11) 



+ 



P 



5G_ _ 6F_ n_G_ 

(5m \ 6p ) (5m y bp 



dx 
dx 



4.3. Derivatives of H. We shall use 
fives. First, 



\(5m' / de 
d 



and compute variational deriva- 
H{m. + e 5m, p) 



e=0 



e=0 



P 



(m + e (5m) dx 



—mdmdx 
P 



Hence 5H/5m = {l/p)m. Next 
d 



6H , 

-T-^^P 
6p 



de 
d_ 
Te 
1 
2 



H{m, p + eSp) 



e=0 



1 I ^-^m^dx+ [{p + e5p)W{p + e5p) 

2 J p + edp J 

^Spm^ dx + j {W{p)+ pW'{p))5pdx 



dx 



Hence 5H/5p = -(l/2)mVp2 + + pW'{p). 

4.4. Computing the Equations of Motion. We start with -F(m, p) 
J mdx. In this case, 6F/5p = and 

SF 



(5m 



(5m 



5m dx 



so 5F/5m = 1. (There must be a more technically accurate way of writing 
this, but for now this will have to do.) We are now in a position to apply 
the general evolution equation 

F{m,p) = {F,H} {m,p). 
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We have 



J dfindx = J m — {l)xdx 
(12) -/ 



m(l)(— ^ dx 

P J X 

2 



1 



p{^)[-i;— + W{p)+pW'{p)] dx 



In the above expression, we have hung on to the symbol (1) everywhere to 
remind ourselves that this is not actually the number one but instead an 
element of the space g*, i.e. a one-form density on M. When applied to 
vector fields, the map (1) acts like the number one. Otherwise we make no 
assumptions regarding its spatial derivatives and charge ahead by integrating 
the first term in ()12|) by parts to obtain 

2^ 



j m^{l),dx = - j {l)dx. 



We get out from underneath the integral sign and also substitute m = pit 

(13) dt {pu) = - {pv?)^ - puux + puua; - px- 

Here p = p'^W'{p) is the pressure. A quick computation reveals that 

Px = P {W + pW')^. 
Thus our final result for the momentum equation is 

(14) dt{pu) + {pu^ +p),^ = 0. 

Next we take F{m,p) = J pdx. In this case, 6F/6m = and just like 
before, 

'6F 



^p^^Pj = j ^Pdx, 

so we shall write 5F/5p = 1, where we mean essentially the same thing by 
"1" as we did earlier. Now the only term from the bracket {F, H} that 
survives is the one involving 5F/5p. Hence 

j dtpdx = j p^^{l)xdx = j m(l)a; dx = - j uix{l)dx. 

Now substituting m = pu and pulling things out of the integral sign we have 



(15) dtp + {pu}^ = 0. 

Equations ((T3|) and (fTl)) give us the system 
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